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The main purpose of the present paper is to analyze the conductance of an A/'-site quantum 
dot with multiple semi-infinite leads. First, we obtain a simple conductance formula that contains 
only the local density of discrete eigenstates and the local density of states of the leads, where the 
discrete eigenstates consist of the bound states, the anti-bound states, the resonant states and the 
anti-resonant states. In other words, the conductance is given by the sum of all the simple poles in 
the complex energy plane. To our knowledge, this is the first time the effect of resonances on the 
conductance is shown exactly. Second, by using the above conductance formula, we show that the 
asymmetry of the Fano conductance peak arises from three origins, namely, from the interference 
between a resonant state and an anti-resonant state, between a resonant state and a bound state, and 
between two different states. Finally, we microscopically derive three types of the Fano parameter 
from the local density of discrete eigenstates, whereas in previous studies in the literature, the 
Fano parameter mainly has been used phenomenologically to describe the shape of an asymmetric 
conductance peak. 

PACS numbers: 72.10.Bg; 73.21. La; 73.22.Dj; 73.23.Ad 



INTRODUCTION 



The electron conduction in nano-scale systems has 
been studied extensively in recent years |l|, 0, 
i, 0, i, i, [Tl, 12, 13, 14, 15, ^,m^^ The 
resonant transport is one of its interesting phenomena, 
where resonant states affect the conductance in its ballis- 
tic transport regime. Resonance is an intrinsic feature of 
open systems (13,1181, 20, 21, 22, 23, 2^,^25", ^,^,^28, 29, 
30, '31', '321,(331, (33, (35!, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 
46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 
62, 63, 64, 65 i6Q,i67|,i6a,l6a,|7Q,i71|,i7a^ 
78, 79, 80, 81]. When we use nano-devices, we inevitably 
attach leads to them. Hence the devices are always open 
systems and have resonant states; an electron comes into 
the device through a lead, is trapped in the confining po- 
tential of the device for a while with a finite lifetime, and 
may come out through another lead. This resonant con- 
duction of mesoscopic systems has been intensely studied 
experimentally H, [TO, [Hi, bJ, [jj, [jj, ll5|, [161, [jj]; for ex- 
ample, the Fano resonance fias attracted mucfi attention. 
The theory of these open systems, however, has not been 
developed so much as applicable to the computation of 
the conductance of nano-scale systems. 

We here consider a class of open quantum-dot sys- 
tems, where all semi- infinite leads are attached to one 
site of a general A^-site quantum dot. For this particu- 
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lar model, we rigorously transform the Landauer formula 
into a simple conductance formula expressed in terms of 
the discrete eigenstates, that is, the bound states, the res- 
onant states, the anti-resonant states and the anti-bound 
states. To our knowledge, this is the first time the effect 
of resonances on the conductance is shown exactly. 

Based on the simple conductance formula, we next dis- 
cuss the symmetry of resonance peaks in the conduc- 
tance. We are particularly interested in an asymmet- 
ric conductance peak, namely the Fano effect [82|- In 
the simplest theory of resonance scattering, a peak ob- 
served in, say, the scattering cross section would have a 
symmetric Breit-Wigner shape of a Lorentzian. In fact, 
some of the peaks are asymmetric. Fano proposed a the- 
ory explaining the asymmetric shape [82]. The asym- 
metric resonance peak has been thereby referred to as 
the Fano resonance. In 2002, K. Kobayashi et al. ob- 
served Fano resonance peaks in the conductance through 
an Aharonov-Bohm system with a quantum dot [§, as 
well as through a T-shaped quantum dot [HI, [lH . Asym- 
metric Fano peaks were clearly observed in the conduc- 
tance. 

It is a conventional understanding that the Fano effect 
arises from the coupling of continuous states in the leads 
and discrete states in the device [qI, [HI, 13 • I^ con- 
trast, we here stress the importance of the interference 
between resonant states [17] as well as between a resonant 
state and a bound state when we consider the Fano con- 
ductance peak. We show that the complex eigenvalues of 
the resonant states of the whole system, the quantum dot 
with the leads, form the asymmetric conductance peak. 

The present paper is organized as follows. In Sec. [Ill 
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we review the theory of resonant states in open quan- 
tum systems. In Sec. IIIH we express, for an N-site open 
quantum-dot model, the retarded and advanced Green's 
functions in terms of the discrete eigenstates. We thereby 
derive a conductance formula consisting only of the local 
density of discrete eigenstates and the local density of 
states of the leads. In Sec. IIVI we show that the asym- 
metry of the Fano conductance peak arises from the in- 
terference between the resonant states as well as between 
a resonant state and a bound state. 



II. RESONANT STATES 

As a preparation for the main part of the present pa- 
per, we review in this section mathematics of the resonant 
state as an eigenfunction of the Schrodinger equation 
It is rather common to define a resonant state as a pole 
of the S matrix. In fact, there are two ways of defining 
the resonant state. The definition based on the S matrix 
may be called the indirect method [83| . We here use the 
direct method of its definition; that is, we describe it as 
an explicit eigenfunction of the Schrodinger equation [H, 

H El, [221, |23|, m i25|, m, [13, m, S, MM. m Hi . 

Suppose that we have a scatterer with several semi- 
infinite leads attached to it. For simplicity and concrete- 
ness, we hereafter restrict ourselves to the tight-binding 
model for the lead Hamiltonians. The total Hamiltonian 
is of the form 



(1) 



where is the one-body Hamiltonian of the scatterer 
(namely, the dot Hamiltonian), is the Hamiltonian of 
a lead a, and is the coupling between the dot and the 
lead a. We assume that 



Ha 



CO 



{\Xa^l){Xa\^\Xa){Xa^l\)- (2) 



Therefore, the energy and the wave number k of in- 
coming and outgoing electrons are related through the 
dispersion relation 



Eh 



-2t cos k. 



(3) 



We can define the resonant state as a solution of the 
Schrodinger equation for the whole Hamiltonian H under 
the boundary conditions that the wave function has only 
out-going waves away from the scatterer [l8|, [lO, [2l|, [221 . 
The condition is often called the Siegert condition [21]. 
More specifically, we seek discrete and generally complex 
eigenvalues En of the whole system 



H\^n)=En\^n). 
{^n\H = En{^nl 



(4) 
(5) 
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FIG. 1: (Color online) (a) Distribution of the eigen- wave- 
numbers kp of the bound states (black crosses), /cp^ of the 
resonant states (blue crosses), k^ of the anti-resonant states 
(green crosses) , and kf" of the anti-bound states (red crosses) 
on the complex wave-number plane, (b) Distribution of the 
eigenvalues Ep of the bound states (black crosses) , Ei^^ of the 
resonant states (blue crosses), E^ of the anti-resonant states 
(green crosses), and Eq^ of the anti-bound states (red crosses) 
on the complex energy plane. The upper and lower halves of 
the k plane respectively correspond to the first and second 
Riemann sheets of the E plane. A branch cut —2t<E<2t 
accompanied by two branch points E = zL2t connect the two 
Riemann sheets. 



in the first Brillouin zone — tt < Re/c < tt under the 
Siegert boundary condition as [3, [85|, [s^ 



{Xali^n) = {i^n\Xa) OC e^^-l^"' 



(6) 



for Xa on any lead a, where is the right- 

eigenfunction and {ijjn \ is the left-eigenfunction [23, [25|, 
[2i,[23,[2i,[2i,[30,[3i|,[32,[83S (Note that (V^^l^ 7^ IV^n) 



3 



in general.) The thus-obtained eigen- wave- number 



as weh as the corresponding eigenenergy 



En 



-2t cos kn 



(7) 



(8) 



are generally complex numbers. Note here that we have 
two Riemann sheets of E for the entire complex plane of 
k (Fig. A branch cut -2t < E < 2t with two branch 
points E = ±2t connects the two Riemann sheets. 

The discrete eigenstates thus obtained are classified as 
follows (Table [J and Fig. [T]). First, the eigenstates with 
/^n > are necessarily on the imaginary axis Re /c = or 
on the edge of the Brillouin zone Rek = n. (In systems 
with continuous space, the bound states exist only on the 
imaginary k axis; the bound states on the line Rek = n 
appear because the leads of the present system are lattice 
systems.) By putting /^^ > in Eq. ([6]), we see that the 
eigenstates are in fact bound states. Hereafter, we use the 
subscript p and the superscript 'b' for the bound states 
as in kp and E^. The bound states with k^p = have 
real negative eigenenergies E^ < —2t while the bound 

states with k^^ = tt have real positive ones E^ > 2t. 

Next, the eigenstates in the fourth quadrant of the k 
plane are referred to as the resonant states. Hereafter, 
we use the subscript / and the superscript 'res' for the 
resonant states as in /cp® and E]^^. The corresponding 
eigenenergies are in the lower half of the second Riemann 
sheet of the E plane: E-^Y^ < 0. 

Third, the eigenstates in the third quadrant of the k 
plane are referred to as the anti-resonant states. (In the 
context of the condensed-matter physics, some refer to 
a resonance in the form of a dip of the conductance as 
an anti-resonance. In the present terminology, this is 
just another resonance, different from the anti-resonant 
state here.) Hereafter, we use the subscript m and the 
superscript 'ar' for the resonant states as in kf^ and E^. 
The corresponding eigenenergies are in the upper half of 
the second Riemann sheet of the E plane: Ei^ > 0. A 
resonant state and an anti-resonant state always appear 
in pair. The states of a pair are related to each other as 



m) = {i'nK and {r^\ = \i;r)^ 



K 



, or 
and 

{Err, or 

E/r and Ei^; 



Z? re 



(9) 
(10) 

(11) 



We refer to a pair of the resonant state and the corre- 
sponding anti-resonant state as a resonant-state pair. 

Some systems have additional states on the negative 
part of the imaginary k axis or on the negative part of 
the edge of the Brillouin zone Rek = tt. Such states 
often appear when resonant and anti-resonant states of a 
pair collide on the axes. We refer to them as anti-bound 
states [89,] and use the subscript q and the superscript 
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FIG. 2: The open quantum N-site dot with multiple leads 
that we consider in present paper. 



'ab' as in k^^ and E^^. Anti-bound states possess real 
eigenenergies but on the second Riemann sheet and still 
have properties of the resonant states such as diverging 
wave functions. 



III. RESONANT SPECTRUM ANALYSIS OF 
AN OPEN QUANTUM iV-SITE DOT 

In the present section, we discuss an A^-site extension 
of the Friedrichs-Fano (Newns- Anderson) model [82|, [H, 
H 112, [H, [H, [95|. We derive a remarkably simple 
conductance formula for the model. The formula con- 
tains only the local density of discrete eigenstates and 
the local density of states of the leads. 

We consider a one-body Hamiltonian of an A'-site dot 
with semi-infinite leads {a} attached to it (Fig. [2]); 



H = Hd-\- ^ {Ha + i^d,a) 



(12) 



with 



N-l 

H^ = ^e,\d,){di\ 

0^i<j^N-l 



\d,){di\).. 



(13) 



Ha 



{\Xa + l){Xa\ + \Xa){Xa + l\), (14) 



Xc,=0 



-ta{\Xa=^){do\ + \do){Xa =0|), 



(15) 



where 6:^, v^j, t and are all real parameters with 
Vij = Vji. The Hamiltonian i^d is the tight-binding 
Hamiltonian of the A'-site dot, while is the tight- 
binding Hamiltonian of the one-dimensional semi-infinite 
lead a and H^ ^ is the hopping between a site on the 
central dot and the end site = ^ oi the lead a. We 
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TABLE I: Classification of tfie discrete eigenstates. 



Bound states 
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first Riemann sheet 
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first Riemann sheet 


E^ > 2t 


Anti-bound states 
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attach all the leads to the single site do of the dot; this 
is due to a technical requirement that appears below. 
The system is an N-site extension of the Friedrichs-Fano 
model [82, sl, H S IH, H [H . We will obtain the 
conductance Q(^^f^{E) from the lead a to the lead (5 in 
the form 



Gmax 



1 ± 



V Pleads (^)y 



(16) 



where 



Peigen 
Pleads (^) 



- E 



27r J-^ E - En 

1 



(17) 



Here Peigen(^) is the local density of discrete eigenstates 
of the whole system on the site do and pieads(^) is the 
local density of states of the lead Hamiltonian "^^Ha- 
The subscripts I and m respectively denote sets 

of the bound states, the anti-bound states, the resonant 
states and the anti-resonant states. 

Let us describe the derivation of Eq. ([16]) hereafter. We 
can obtain the exact expression of the scattering states of 
the system (p!2]) . namely the Friedrichs solution [90] IV^^), 
with the eigenvalue Ek = —2t cos k; see Appendix[Al The 
completeness with respect to the scattering states is given 
by m 
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(18) 



where is a bound state and is a scattering state 
given in Appendix [A] and we used the notation 



FIG. 3: (Color online) (a) The integration contour C^z (pur- 
ple curve) for the retarded Green 's function G^(E) and (b) 
the integration contour Cbz (cyan curve) for the advanced 
Green's function G^{E), with the circular (gray) contours 
extracting bound states in the complex wave-number plane. 



(19) 



We first express the retarded and advanced Green's func- 
tions in the spectral representation; 



where the integration contours (7^2 and (7^2 cover the 
Brillouin zone as indicated in Fig. [3l 



p 

G\E) = Y^ 



E 



E^ 



E 



E^ 



CA 27r E-Ek 



Next, we replace the integration contours C^^ ^^'^ ^bz 

(20) 

with the ones shown in Fig.[4l Then the retarded Green's 
function G^{E) acquires the residual integrals of the res- 
onant states /cp®, which lie in the fourth quadrant, while 
(21) the advanced Green's function G^{E) acquires the resid- 
ual integrals of the anti-resonant states /c^, which lie in 
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FIG. 4: (Color online) (a) The integration contours C^{i<io) 
and Cf{K,o) (purple curves) for the retarded Green's func- 
tion G^{E), modified for extracting the resonant states (blue 
crosses) in the complex wave-number plane, (b) The integra- 
tion contours C±{k,o) and C^{k,o) (cyan curves) for the ad- 
vanced Green's function G'^(E), modified for extracting the 
anti- resonant states (green crosses). 



the third quadrant: 



We then have 



G^(^) 



E 



E^ 



lim / 

^+°° Jc' 



dk 



b\/J,b| ^ 



27r E 



(24) 



E^ 



E 



lim / 



CA(/^o)+CA(/^o) 



27T E-Ek 



(25) 



Here C^{hZo) indicates the sum of the paths parallel to the 

real axis and C^{hZo) the sum of the paths perpendicular 
to the real axis including the contributions from the anti- 
bound states. Note that hzo of the modified integration 
contour must be positive and greater than the imaginary 
parts of all the resonant eigen- wave- numbers. 

At this point, we sum up the retarded and advanced 
Green's functions; 



G^{E)^G^{E) 



=2E 



E-E^ 



E 



res I 



E 



^res 



E 



E-E?- 



Ko^oo J^R(^^)_^(-,A(^^) 2tt E-Ek 



lim 



Cf{Ko)+Cf{Ko) 



27r E — Ek 



(26) 



The sum of the contributions of the integration contour 
C^^kq) and C^^kq) is equal to the contribution of the 
bound states and anti-bound states except for the sign; 



lim 



dk \^l){^l\ 
27r E — Ek 



E 



iv^f)(^f I 



E 



E^ 
j^p 



E 



(27) 



On the other hand, we proved that the contributions of 
the parallel integration contours C^{nQ) and C^{nQ) van- 
ish for the states on the central dot; i.e. for \di) and \dj) 
with any i and j, we have 



lim 



Cf{Ko) + Cf{Ko) 



27r E-Ei, 



0. (28) 



27r E-Ek 



res I 



E - Ef 



JC{k=k^) 



2n E-Ek 



E-E- 



(22) 
(23) 



See Appendix [B] for the proof. Equation ([28|) does not 
seem to hold if the semi-infinite leads are not attached to 
a single site of the central dot. This is why we focused 
on the present system (p!2]) . 

Thus we find that the sum of the retarded and ad- 
vanced Green's functions is equal to the contributions of 
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FIG. 5: (Color online) Summation of the modified integration 
contours for the retarded and advanced Green's functions in 
the complex wave-number plane results in the sum of the 
contributions of all discrete eigenstates containing the bound 
states (black crosses), the resonant states (blue crosses), the 
anti-resonant states (green crosses) and the anti-bound states 
(red crosses). 



only the discrete eigenstates for the states on the central 
dot, {\d,)}, (Fig. El); 



{di\ {G^{E) + G^{E)) \dj) = {di\A{E)\dj), (29) obtain the final solution as 



see Appendix [Cl for derivations of the self-energies. (The 
relation ([3T]) holds for arbitrary sites, not restricted to the 
states of the central dot, {\di)}.) Equation ([29]) shows 
that the real part of the Green's function is given by 
the discrete eigenstates, while Eq. (|3T]) shows that the 
imaginary part of the Green's function is given by the 
inverse of the van Hove singularities at the branch points 
E = ±2t 0. 

The simultaneous matrix equations ([29]) and (|3T]) re- 
sults in the matrix Riccati equations 

{d,\{G^{E)mE))G^{E) 

+G^{E) [2 + {ir{E)) A{E)] + A{E)} \dj) = 0, (36) 

{d,\{G^{E){-iTiE))G^iE) 

- [2 + A{E) {-iT{E))] G^{E) + A{E)} \dj) = 0. 

(37) 

The solution gives each Green's function in terms of the 
contribution of the discrete eigenstates, A(£^), and the 
contribution of the branch-point singularities, r(^). Us- 
ing the fact that {do\r\do) is the only non-zero element of 
the matrix {di\T\dj) for the present system (p!2l we first 
solve the above equations for i = j = 0, then for z = 
with general j and for j = with general z, and finally for 
general i and j. Denoting {di\G^/^{E)\dj) = G^J^{E), 
{di\A{E)\dj) = A,j{E) and {d,\r{E)\dj) = r,,(E), we 



where 



A(^)- E 



n^p,q,l,m 



E — En 



(30) 



We also use the fact that the difference between the re- 
tarded and advanced Green's functions is generally given 

by i 



Gl{E) 



where 



Aij{E) , Aio{E)Aoj{E)l±^l-n{EY 



Aoo(i?) 



2in 



(38) 



Aij{E) Aio{E)Aoj{E)l±^l-n{Ef 



G^{E) - G^{E) = iG^{E)T{E)G^{E), (31) 



Aoo(^) 



n{E) = Itoo{E)Aoo{E). 



2in 



(39) 



(40) 



where 



T{E)^i{V,^{E)-V,%{E)) Mo)(do| 

= ^r(")(£) 

a 

with the self-energies of the semi-infinite leads 
'taV E -iV4.t^ - E^ 



rv ^ ^ 



(32) 

(33) 
(34) 



The sign in front of the square root of Eqs. (|38|) and (|39|) 
is chosen according to the rule given in Appendix [Dl 

Using the Fisher-Lee relation [97], we arrive at the con- 
ductance Qa^f^{E) from the lead a to the lead /3 in the 
form 

2e2 



Ga^piE) = ^Tr (t(-^\E)G^{E)T^'''>{E)G^{E)') 

= ^r^^,\E)GUE)r^^HE)GUE) 

Gmax 
a^f3 



1 ± - n{E)^ 



and 



where 



r(«)(^) ^ ( ) ^4fi-E^\do){do\; (35) 



Gmax 



2e^ I 2tatft 



(41) 



(42) 
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FIG. 7: A point contact do with two leads. 



conductance is exactly given in terms of the sum of sim- 
ple poles of the discrete eigenstates. 



FIG. 6: (Color online) The local density of the discrete eigen- 
states, peigen{E) (bluc curvcs) , and the local density of the 
lead states, pieads{E) (red curve). The plots are for a two-site 
dot with two leads. The parameter values are ti/t — t2/t — 1, 
Eo/t — 5, Si/t — and VQi/t — vio/t — 1. The gray verti- 
cal lines indicate the two bound states (^i = —5.20735 . . ., 
E2 — 5.55278 . . .), as well as the real part of the resonant-state 
pair (^r^^/^r = -0.172712 . ..^^0.0689508 . . .). The structure 
of Peigen {E) around the resonant-state pair is continuous for 
real E. The shaded areas indicate the outside of the energy 
band, |^| > 2t. 



is the maximum possible conductance from the lead a 
to the lead (3. (In the transformation of Eq. (|4T]) . we 
again used the fact that the matrix F^"^ has only the 
(0,0) element for the present system (p!2]) : see Eq. (|32]) .) 
Equation (|^ gives a remarkably simple formula 



Gmax 



1 ± 



\p\eads{E)J 



(43) 



where 

Peigen 
Pleads 



Aoo(^) ^ i 
27r 



27r 
1 



- y 

n^p,q,l,m 



{do\^l^n){'lpn\do) 
E — En 



Trroo {E) Y.a (ta /tf TtVW^E^ ' 



(44) 



Here Peigen(^) is the local density of discrete eigenstates 
of the whole system H on the site do, whereas pieads(^) 
is the local density of states of the lead Hamiltonians 
i^Q,, which has the van Hove singularities at the band 
edges E = ±2t. Note that Peigen(^) has singularities at 
the discrete eigenvalues, whereas pieadsiE) has singular- 
ities at the branch points. The conductance itself has 
singularities due to the discrete eigenstates but not due 
to branch points. We exemplify peigen{E) and pieads(^) 
in Fig. [6] for a two-site dot with two leads with h/t = 
t2/t — 1, eo/t = 5, ei/t — 0.5 and voi/t = vio/t = 0.5. 

To summarize the present section, we reveal the effect 
of resonances on the conductance explicitly and rigor- 
ously. To our knowledge, this is for the first time the 



IV. QUANTUM INTERFERENCE EFFECT OF 
DISCRETE EIGENSTATES 



In the present section, we argue that the Fano conduc- 
tance arises as a result of interference between discrete 
eigenstates. The conductance formula (|^3|) has a square 
of the local density of discrete eigenstates. Therefore, we 
have crossing terms within a resonant-state pair (between 
a resonant state and an anti-resonant state), between two 
resonant-state pairs (two sets of a resonant state and an 
anti-resonant state), and between a resonant-state pair 
and a bound state. We show in the present section that 
discrete eigenvalues decide the symmetry or the asymme- 
try of the conductance peaks in addition to the location 
of the conductance peaks, using several examples. We 
thereby derive the Fano parameter microscopically. In 
Subsecs. A, B, and C of the present section, we consider 
the system ([12]) with the following restrictions: only two 
leads a = 1,2; the coupling ti = ^2 = t; the number of 
sites in the dot N = 1,2,3. We consider the effect of 
changing in Sec. IIVDI Throughout the present sec- 
tion, we computed the conductance using the Fisher-Lee 
relation (jH]) and obtained all discrete eigenvalues solving 
Eq. ([C30l) . 



A. Point contact system: N = 1 

First we show the conductance as well as the discrete 
eigenvalues of the one-site dot with two leads, namely the 
point contact shown in Fig. [71 There are only two bound 
states and no resonant state. We plot in Fig. [8] the con- 
ductance with the eigenvalues of the two bound states for 
eo/t = 0, 1, 1.5, 2, 2.5. The conductance of the point con- 
tact has no peculiar behavior such as the Breit-Wigner 
peak or the Fano peak. Upon increasing the potential 
£0, the eigenvalues of the two bound states move away 
from the branch points E = ±2t. This decreases the con- 
tribution of the local density of the discrete eigenstates 
Peigen(^) and hcucc deflates the conductance gradually. 
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FIG. 8: The energy dependence of the conductance (the left 
axis) and the discrete eigenvalues of the bound states (the 
right axis) for the one-site dot with eq/I = 0, 1, 1.5, 2, 2.5. 
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FIG. 9: A two-site quantum dot with two leads. 



B. T-shaped quantum-dot system: N — 2 



We next show the conductance and the discrete eigen- 
values of the two-site quantum dot with two leads, 
namely a T-shaped quantum dot shown in Fig. [9l This 
system is a minimal model that possesses a resonant- 
state pair (a resonant state and the corresponding anti- 
resonant state) and may be directly related to Fano's 
original argument [82]. We plot in Fig. [10] the conduc- 
tance, the eigenvalues of the two bound states, E\ and 
£^2 : the eigenvalues of the resonant-state pair, E"^^^ 
and for £0/^ = 0, 1, 3, 5, £1 = and voi/t = vio/t = 
1. 

We have a Breit-Wigner dip for sq = 0^ but for 7^ 
we have an asymmetric peak, namely the Fano conduc- 
tance peak. Maruyama et al [98] claimed that the asym- 
metry of the conductance peak of the T-shaped quantum 
dot is proportional to Sq. We here discuss the asymme- 
try from the viewpoint of interference among the discrete 
eigenstates. 

The conductance formula (|^3|) contains the square of 
the sum over the discrete eigenvalues of the form 



n{Ey 



f PeigenjE) 
V Pleads (£) 



^pair 



iE)y 



(Pleads (-E))' 



(45) 



0.8 
0.6 



a 0.4 



o 0.2 











X 




X 


K 






>; 













(c) 



1 

:^ 0.8 

a 0.4 
-a 
o 

o 0.2 



-2 2 

Elt.EJt 



-2 2 
Elt.EJt 









8o/^ = 5 









¥r- 










1 









0.6 



0.2 
' 
-0.2 



-0.6 



0.6 



0.2 



-0.6 



FIG. 10: (a) The £0 dependence of the conductance (the left 
axis) and the discrete eigenvalues (the right axis) for the two- 
site dot with (a) so/t = 0, (b) so/t = 1, (c) so/t = 3 and (d) 
So/t = 5. Here we fixed Si/t = and v^i/t — vio/t — 1. 
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FIG. 11: (Color online) The quantities Vl(EY (gray curve), 
VH'iEf (chained green curve), Qp^^'(^)^ (broken red curve 
and r2^-P^^^(^) (dotted blue curve), defined in Eqs. (|45 l) - (|5T]) . 
plotted with the conductance (solid black curve), Eq. (|4T]) . 
or Eq. (|43p . The system is the two-site quantum dot. We 
fixed £o/t = 5, exji — and v^xji — vio/t = 1. The gray 
vertical line indicates the real part of the resonant eigenvalue 
El^"" = -0.172712.... 



where 



p\E) 



^ 1 {do\^^){^^\do) 



(46) 



2n E- E'' 



27r E- E'' 



(47) 



^^{Ef 



(Pleads (-E)) 
(-pPair(-^))2 

(pieads(-E^))^ 



2 ' 



^b-pa.(^) ^2p^i^)pPai-(i^) 



(Pleads(-B)) 



(49) 

(50) 
(51) 



The second quantity (|5Q|) contains a crossing term be- 
tween the resonant state and the anti-resonant state. The 
third quantity ([5T]) contains crossing terms between the 
resonant state and a bound state as well as crossing terms 
between the anti-resonant state and a bound state. We 
can see in Fig. [11] that the asymmetry of the conduc- 
tance peak comes partly from the asymmetry of the term 
^pairj^^j and partly from the crossing term 1]^"P^''"(£;). 
The quantity r!!^{E) is almost symmetric. 

In order to derive the Fano parameters for the asymme- 
try of the two terms ^^^^^^{E) and ^^'"^^^^{E) microscop- 
ically, we expand the terms ([5Q|) and ([5T|) in the neigh- 
borhood oi E = El^^ = by using the normalized 
energy 



E 



E - E'/ 
~\Ef^ 



We first rewrite p'^^^^{E) in the forms 



p^^''{E) 



E - (£;re« + iEl""^) 



+ c.c, 



(52) 



(53) 



where we express the coefficient of the local density of 
the resonant state with the amplitude N and the phase 



_ {do\r''){r''\do) 



(54) 



Since the conductance formula (|43p is given in the form 



Goc 



n{Ef 



1 T VI - ^{Ef 



(48) 



the symmetry or the asymmetry of the quantity Vt[EY is 
directly reflected on the symmetry or the asymmetry of 
the conductance peak. Equation (|45|) therefore implies 
that the symmetry or the asymmetry of the conductance 
peak is strongly affected by crossing terms, or the in- 
terference between states with discrete eigenvalues. We 
hereafter show that the Fano conductance peak arises 
from two types of the interference, or two types of cross- 
ing terms. First, we have a crossing term within the 
resonant-state pair, or the interference between the reso- 
nant state and the anti-resonant state. Second, we have a 
crossing term between the bound states and the resonant- 
state pair. 

We compare in Fig. [11] the following quantities: 



Note that this is generally a complex number because the 
left-eigenvector (?/;^®^| is not generally Hermitian conju- 
gate to the right-eigenvector |?/^^®®) for a resonant state 
(see Eq. (j9j)). We then rewrite the local density of the 
resonant-state pair in the form 



,pair(^) ^ ^, {E-Er)cose^\Er\smt 



{E - El'^^Y ^ 
N sin 6> + ^ cos i 



\El 



or 



1+^2 



i + £;2 



where 



tan^. 



(55) 



(56) 



(57) 



The parameter (|57|) controls the asymmetry of the 
term (|5Q]) and hence may be called the Fano parameter. 
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although Eq. ([56l) is diflFerent from the form originahy 
derived by Fano [82,]: 



G 



1 + ^2 



(58) 



The asymmetry caused by the above interference between 
a resonant state and the corresponding anti-resonant 
state may be missing from Fano's argument. 

On the other hand, the crossing term ([5T]) produces 
asymmetry of Fano's original form ([58|) . In order so 
see this, we approximate the local density of two bound 
states as 



A) ( T^re 
P {^1 



(59) 



in the neighborhood of £^ = E^^^. We therefore have the 
crossing term between the resonant-state pair and the 
two bound states as 



where 



- sE + tE'^ 

1 + ^2 ' 



_ p'^iEr) 



(60) 

(61) 

(62) 
(63) 

In order to derive a Fano parameter q^-P^^^ that controls 
the asymmetry of the term Q^~^^^^{E), we extract the 
form on the right-hand side of Eq. ([58|) by putting 

r^sE^ tE'^ 

ihTP 



^=^^cos^ + /(^r)sin^, 
t = p^\eI^') cos 0. 



/"^b-pair ^^A^ 

a^b^ . — ^, (64) 



We obtain the Fano parameter gt)-pair solving the equa- 
tion 



5 = 



(65) 



and choose the solution with the same sign as s. This 
controls the asymmetry of the term (|5T]) . a Fano param- 
eter that is different from the one given by Eq. ([57|) . but 
that conforms to Fano's original form ([58|) . 

We show in Fig. [12] how the two Fano parameters q"^^^^ 
and gr^-pair (^gpg]2(^ the system parameter ^o- In the 
particular case of Fig. [121 q^'^^^"^ tends to dominate over 
^pair increase the system parameter ^o- This is 

in coordination with the decrease of We can see 

in Eq. (|62|) that a small imaginary part \El^^\ causes a 
particularly strong asymmetry of the term Q^'^^^^^E). 
This is indeed demonstrated in Fig. [TOl where, as we in- 
crease the asymmetry rapidly develops while the the 
resonant eigenvalue approaches the real axis. Inciden- 
tally, the present system has the particle-hole symmetry 
E ^ -E for So = ei = 0, and hence gP^'"" gb-pair ^ 
for which the resonance peak takes the form of a sym- 
metric Lorentzian as shown in Fig. [TOlfa). 




FIG. 12: (Color online) The Fano parameters q^''''' (blue 
curve) and gr^-P^^^ (red curve) for the two-site dot, plotted 
with the imaginary part of the resonant eigenvalue El^^ . Use 
the right axis for the Fano parameters and the left axis for 
the resonant eigenvalue. (Note that the horizontal axis does 
not start from sq = 0. We omitted the part < £o < 2 
to avoid confusion because the structure of the spectrum is 
drastically different in that region.) We fixed Si/t = and 
voi/t — Vio/t — 1. 



Ed 




FIG. 13: The three-site quantum dot with two leads. 



C. Three-site quantum- dot system: N = 3 

Third, we discuss the conductance of the three-site 
quantum dot with two leads shown in Fig. [131 This sys- 
tem have two resonant states for some parameter val- 
ues. This situation was not considered in Fano's argu- 
ment We show in Fig. [TH the conductance, the 
eigenvalues of the two bound states, E^ and ^2 ^ well as 
the eigenvalues of the two resonant-state pairs, El^^, Ef^^ 
Ef^ and Ef, for si/t = -1.5, -1, -0.5, with eo/t = 0, 
e2/t = 0.5, voi/t = vio/t = 0.8, vo2/t = V2o/t 0.5 and 
vi2/t = V2i/t = 0.4. Upon increasing the parameter si, 
the conductance dip that is generated by the resonant 
state on the left-hand side, El^^^ approaches to the other 
conductance dip that is generated by the resonant state 
on the right-hand side, £^2^^. Then the latter conductance 
peak develops strong asymmetry. 

For the present system, we have yet another Fano pa- 
rameter due to a crossing term between one resonant- 
state pair and the other resonant-state pair. The con- 
ductance formula (|^ contains the square of the sum 
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FIG. 14: The conductance (curve for the left axis) for the 
three-site dot with (a) Si/t = —1.5, (b) £i/t = —1.0, (c) 
exji — —0.5 and (d) exji — 0, plotted with all the discrete 
eigenvalues (crosses for the right axis). We fixed £o/t = 0, 
s^jt — 0.5, v^xji — vio/t — 0.8, VQ2/t — V2o/t — 0.5 and 

Vl2/t — V2l/t — 0.4. 



over the discrete eigenvalues of the form 

[p\E)+pr{E)+pr{E))' 



VL{Ef 
where 

p\E) ^ 



E 



(Pleads 
1 (rfo|^b)(^^|^^) 



27r E- 



(66) 



(67) 



p=l,2 

. 1 {d^wr){^r\do) 



1 (doi^r)(^rido) 



'27r 



E- 



^res 



E-Ef 



27r 

for / = 1,2. 

We compare in Fig. [15] the following quantities: 

\2 



\Ef 
\E) 

l]p^^^-p^^^(£;) 



^b-pair^ 



(pleads(-B))' 
(Pleads(-E))' 

2p^{E)pr^{E) 

(Pleads (-B))^ 

2pr\E)priE) 



for I = 1,2, 



for Z = 1, 2, 



(Pleads (-B))' 



(68) 

(69) 

(70) 
(71) 
(72) 



We can see that the following three terms are asymmet- 
ric: first, 1^2^^^ which contains the crossing term be- 
tween the resonant eigenstate ip2^^ and the anti-resonant 
eigenstate ^/^f^; second, 0.2 '^^^^ i^)^ which is the cross- 
ing term between the bound states (7/^5^,7/^2) 
resonant-state pair (7/^2^®,?/^!^); third, (^p^^^"P^^^(E), which 
is the crossing term between the two resonant-state pairs 

(V'r^V'f) and (V'r^V'f)- 

In order to derive the Fano parameters for the asym- 
metry of the three terms, we expand the terms ([7Q]) - ([72|) 
in the neighborhood of £^ = £^^1^ by using the normalized 
energy 



E = 



E 



Tpres 



\E!r\ 



(73) 



We can analyze the terms 0.2^" (E) and q!^''^^"{E) in the 
same way as in the previous subsection. We again use 
the expression 



(doi^r)(^rMo) 



(74) 



Then the Fano parameter controlling the asymmetry of 
the term 1^2^^^ (£^) is given by 



pair 
92 



tan^. 



(75) 



Following the same logic as in Eqs. ([52 ]) -(|65 |) , we obtain 
the Fano parameter that controls the asymmetry of the 
term Q2~^^^^{E) by solving 



b-pair 



0, 



(76) 
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FIG. 16: (Color online) The Fano parameters ^2^^^ (blue 
curve), ^2"^^^^ (^^^ curve) and ^p^^^-p^^^ (purple curve), plotted 
with the difference of the real parts of the two resonant eigen- 
values, — ^ri^- Use the right axis for the Fano parame- 
ters and the left axis for the eigenvalue difference. We fixed 

£o/t — 0, S2/t — 0.5, Voi/t — Vio/t — 0.8, V02/t — V2o/t — 0.5 

and vi2/t = V2i/t = 0.4. 



use the expansion 



(80) 



We then approximately have the crossing term between 
the two resonant-state pairs as 



FIG. 15: (Color online) The quantities Q(E)'^ (gray curve), 
Q'^{E)'^ (green curve), Qf^^'(£;)2 (broken and chained red 
curves), Q^~'^^^'^{E) (broken and chained blue curves) and 
^pair-pair (dottcd purplc curvc) , defined in Eq. (|66l)-(|72l). 
plotted with the conductance (solid curve), Eq. (|4T]) . or 
Eq. (|43p . The vertical gray lines indicate the real parts 
of the resonant eigenvalues E = = —0.211544... and 

E = Elf = 0.721170 . . .. (b) shows the part of (a) around 
E — El2^ with the plots of the conductance (solid curve), 
^{Ef (gray curve), ''(£;) ^ (chained red curve), Q:^{^'''\E) 
(chained blue curve) and Qp^^''"P^^'"(^) (dotted purple curve). 
The system is the three-site dot. We fixed so/t = 0, Si/t = 0, 
S2/t — 0.5, VQi/t — vio/t — 0.8, VQ2/t — V2Q/t — 0.5 and 

Vl2/t — V2l/t — 0.4. 



l]P^i^-P^i^(^) 



with 



2pf"(^)pf"(^) ^ r' + s'E + t'E'' 

{plea.ds{E)f 



1 + ^2 



■ sm ( 



pair- / 771: 
Pi (A 



cos6> + pf^' (^rT) sin i9, 



res\ 
r2 J 



(81) 

(82) 

(83) 
(84) 



We thus have yet another Fano parameter gfPair-pair 
solution of 



where 



^i2 I 
^o/Z7'res^ 



' \ET. 
t = p^'{El^'')cos( 



(77) 

(78) 
(79) 



Next, in order to discuss the quantity 1]p^^^"P^^^(£^), we 



S' (gPair-pair>j- _ _ ^/)^pair-pair _ ^/ = Q. (85) 

We show in Fig. [16] how the three Fano parameters 
^2^^^, ^2"^^^^ and gpair-pair (^gpgj^^^ ^hc systcm param- 
eter El. In the particular case of Fig. [161 the third Fano 
parameter gP^ir-pair greatest in most of the range. 

This may be due to the following reason. The first term 
of for the parameter gP^ir-pair contains the Lorentzian 



pair/T7ires\ / 771 res Tpres\z , zrires 



(86) 
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Therefore, s' grows fast as the resonant-state pair 
approaches the resonant-state pair £^^"1® up until 
This is in contrast to the first 



^res 

i4t 



term of s for the parameter q}^''^^" , which contains 



(87) 



for p = 1,2. This is indeed demonstrated in Fig. [HJ 
where, as we increase £1, the asymmetry rapidly devel- 
ops while the resonant-state pair {El^^^Ef^) approaches 




D. The effect of the hopping energy to 
central dot and the leads 



between the 



Finally, we briefly show the effect of the hopping en- 
ergy ta between the central dot and the lead a. We here 
use the case of the three-site dot with two leads with 
h = t2 ^ t, eo/t 0, ei/t 0, £2/^ 0-5, voi/t 
vio/t 0.8, vo2/t V2o/t 0.5 and Vi2/t V2i/t 0.4. 
For ti = t2 < t/V^, there are three resonant-state pairs 
and no bound states. We have corresponding three sharp 
peaks in the weakly coupled case ti/t = t2/t = 0.1 
as in Fig. \17\ (a). Upon increasing the hopping energy 
ti = ^2, the second peak corresponding to the resonant- 
state pair with the least modulus of the imaginary part 
develops asymmetry. At ti/t = t2/t = 1/a/2, the reso- 
nant and anti-resonant states of a resonant-state pair col- 
lide and become two anti-bound states, which leaves two 
resonant-state pairs. For ti/t = t2/t > 1/a/2, the second 
peak continuously develop the asymmetry. (The anti- 
bound states become bound states before ti = t2 = t.) 



V. CONCLUSION 

We carried out the spectrum analysis of the open quan- 
tum N-site dot with multiple leads. We obtained the 
simple conductance formula (|^3|) in terms of the local 
density of discrete eigenstates (the bound states, the res- 
onant states, the anti-resonant states and the anti-bound 
states), peigen(^), and the local density of states of the 
leads, Pleads (^)- To our knowledge, this is the first time 
the conductance is exactly give by the sum of all the 
simple poles. 

We then showed that the Fano conductance arises from 
the crossing terms of three origins; first between a pair of 
a resonant state and an anti-resonant state, second be- 
tween a resonant-state pair and a bound state, and finally 
between two resonant-state pairs. We also presented mi- 
croscopic derivation of the Fano parameter. 

The analysis in the present paper is applicable only to 
non-interacting systems. It is an interesting and chal- 
lenging problem to generalize the present approach to 
interacting systems. The Kondo effect, for example, has 
been observed in recent experiments on quantum dots 
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FIG. 17: The conductance (curve for the left axis) for the 
three-site dot with (a) ti/t = t2/t = 0.1, (b) ti/t = t2/t = 0.3, 
(c) ti/t = t2/t = 0.6 and (d) ti/t = t2/t = 0.8, plotted with 
all the discrete eigenvalues (crosses for the right axis) The 
gray curves and the gray crosses indicate the conductance 
and the discrete eigenvalues for ti/t = t2/t — 1, the same 
data as plotted in Fig. [M] . We fixed so/t — 0, Si/t — 0, 
S2/t — 0.5, VQi/t — vio/t — 0.8, VQ2/t — V2Q/t — 0.5 and 

Vi2/t — V2l/t - 0.4. 
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and attracts much theoretical interest. The present ap- 
proach may be particularly useful in analyzing the inter- 
play between the Fano resonance and Kondo resonance. 
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APPENDIX A: FRIEDRICHS SOLUTION OF 
THE SYSTEM ([T2]) 

In the present appendix, we solve the Lippmann- 
Schwinger equation for the present system (p!2|) to obtain 
the Friedrichs solution [93] of the scattering states. The 
Lippmann-Schwinger equation may be written down as 



we then have 



\i^lo) = \k,a) + 

where 

Hq =Hd + ^ Ha 

a 

N-1 

= Y.e\d,){di\ 

i=0 



1 



Hi\^la). (Al) 



0<i<j<Ar-l 

^ / '^'^^ ' 



(A2) 



a 

= -E*"/ ^{\k,a){do\ + \<k){k,a\), (A3) 

the state a) is an eigenstate of Hq (more specifically, of 
Ha) with the eigenvalue Ek = — 2tcosA:, and 5 is a posi- 
tive infinitesimal ensuring that the solution is an outgoing 
wave. 

The formal solution of the Lippmann-Schwinger equa- 
tion ()Aip is given in the form 



\4'la)=\k,a) + 



-M,\k,a) 



Ek-H + i5 
k,a)- ^ , .J do). (A4) 



Ek-H + i5 
Using the resolution of unity 

dq 



1 = ^ \d,) {d^ + Y^f ^ \q,P) {q,P\ , (A5) 
i=o p J-^ ^'^ 



\€.a)=%a)-ta ( ^Gfo(^.)|^. 

1 



27r 



(A6) 



where 



1 



Ek-H + i5 



Mi) • (A7) 



In order to transform the final term on the right-hand 
side of Eq. ()A6p . we calculate the following: 



1 



Ek-H-iS 



77 



^ -Ml]- 1 -\q,l3) 



Ek-H -i5 J Ek- Ho- iS 



1 



Ek-Eq-iSy^"' Ek-H-iS 
We thereby have 

1 



\do) ■ (A8) 



{q,f3\ 



Ek-H + i5 
We therefore arrive at 

/N-l 



- \do) = -tpG%{Ek). (A9) 



=0 



" dq tp\q,l3) 
2tt Ek-Eq + iS 



. (AlO) 



We describe in Appendix \C\ how we can calculate the 
Green's function Gf",. 



APPENDIX B: PROOF OF EQ. 



In the present Appendix, we prove Eq. (|28|) . Using the 
expression (jAlOp of the scattering state, we have 

a 

= Y.tjGl{Ek)G^j{Ek). (Bl) 
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We therefore have 



dk {d^\^l){^Pl\dj) 
27r E-Ek 

dk 1 



a 

X (d 



27r£; + t(e^^ + e"^^) 
1 



X (c^ol 



E-H^-^{tJ/t)e'^\do){do 



E-H^-^{tJ/t)e-'^\do){do 



■\do) 



■\dj). (B2) 



where we used Eq. ()C14p for the Green's functions with 
the expression (|C33p for the effective potential. 

On the paths C^{hZo) and C^{hZo)^ we let k = ^in^ 
and integrate with respect to k^. For k = k^ -\- ihco^ the 
element e~*^ grows in the limit a>:o ^ oo in the denomi- 
nators of two of the three factors on the right-hand side 
of Eq. (|B2p . For k = k^ — in^^ the element e*^ grows in 
the limit /^o ^ oo again in the denominators of two of 
the three factors. Therefore the integral (|B2p vanishes 
on the paths C^{i<iQ) and Cy^(/^o) in the limit oo. 
Thus Eq. ([28]) is proved for the system (fT2|) . 



APPENDIX C: THE GREEN'S FUNCTION IN 
THE CENTRAL DOT AND CALCULATION OF 
THE RESONANCES 

In this appendix, we describe the calculation of the 
Green's function G^j{E) for the states in the central dot, 
The calculation utilizes the self-ener gy o f the 
semi-infinite leads [2, [H, Hi, H [ToO, EoH, |io2, lioi, Eol, 
[l05, 106, 107, 108|. Using the expression of the Green's 
function, we also give an equation that gives the resonant 
states. 

The basic statement is the fact 



GfAE) = {di 



1 



(Cl) 

where the thus-defined effective Hamiltonian has de- 
grees of freedom only on the central dot. Below, we will 
review the derivation of the following form: 



Hfs{E) = Ha + V,^{E)\do){do\ 



(C2) 



where 



AT-l 



H^ = Y^ ei\di){di 



=0 

E 

0<i<j<N-l 
2 



vi,m{dj\ + \d,){d.\), (C3) 



(C4) 



Therefore, we can calculate the Green's function by 
inverting an x matrix (|C2p . 

There are several ways of deriving Eq. (|Cip . One way 
is to use the resolvent expansion 

1 _ 1 

E-H^iS ~ E-Ho^iS 

I 1 rr 1 

E-Ho^iS ^E-Ho^iS 

1 .i^.^^H. ' 



E-Ho^iS "E-Ho^iS " E - Ho ^ iS 

where 

Ho =i^d + ^ Ha 

a 

N-1 

= J2£\di){di\ 

i=0 

- E ^iji\di){dj\ + \dj){di\) 

0<i<j<N-l 



(C5) 



- *E E + + \^o.){Xa + 1|) (C6) 

Hi = ^ Hd^a 

a 

= -J2ta.{\Xa=0){do\ + \do){Xa=0\) (C7) 

a 

In calculating (E) defined in Eq. ()Cip . we should note 
the following. Let Hd denote the Hilbert space spanned 
by the states on the central dot, and TYiead denote 

the Hilbert space spanned by the states on the leads, 
{|xck)}. Then we have 



1 



■■\d^) = 



1 



E — Ho -\- 10 E — Hd + to 

i I -x l^^) = i , -A ^a) eHlead, (C9) 

E — Ho ^ 10 E — Ha + ^0 

Hi\di) = -SoiY,ta\Xa = 0) GHlead, (ClO) 

a 

Hi\xa) = -Sa^Mdo) eHd. (Cll) 

That is, the operator {E — Ho -\- iS)~^ , when applied to 
a state either in Hd or TYiead, does not change its Hilbert 
space, whereas the operator Hi switches it. Therefore, 
all terms of odd orders of Hi in the resolvent expansion 
of vanish. All terms of even orders of Hi (except the 
zeroth order) have powers of the following factor: 



Y,{do\Hi\xa=0) 



X {Xa 0|- 



1 



'E-Ho^iS 

X {Xa = 0|i^i|do) 



\Xa = 0) 
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We will show below that this quantity is equal to V^{E) 
defined in Eq. ()C4p . We therefore have 

GfAE) 



{di\ 



{d 



1 



'E-HA+i5 



|do>Kff(^^o| 



1 



which can be summarized as 



E-HA-\do)V^{do\+i5 



E-Hd + i5'^" 

(C13) 



Id,-). (C14) 



The remaining task is to calculate {x^ = 0|(£^ — Hq 
i6)-^\xa = 0) in Eq. (p^ . or 



GLd(^;0) = (x = 0| 



1 



|x = 0), (C15) 



where we fixed the sign in front of the square root so that 
the imaginary part may be negative. Thus the quan- 
tity ()C12p was indeed shown to be equal to V^{E) de- 
fined in Eq. ([U4l) . 

To summarize the above, the retarded Green's function 
given by ()C14p is expressed in the form (|Cip with the ef- 
fective potential V^^ defined in Eq. (|C4p . The Green's 
functions that are used in the expression of the scatter- 
ing state (|AlQp are therefore obtained by inverting the 
N X N matrix {di\{E - H^^{E))\dj) for a fixed value of 
E. Incidentally, the infinitesimal -\-i6 in the denominator 
of the definition (|Cip is not necessary anymore because 
already has an imaginary part. In fact, the advanced 
Green's function is given by fiipping the sign of the imag- 
inary part; 



{di\ 



1 



E-H -i5 



E 



\d,) 
(C22) 



with 



where 



=X 



We then use the resolvent expansion 
1 1 



E-Flead(0)+i(5 E-fflead(l) + i<^ 
1 



E-fflead(l) + i<^ 

x(-t)(|l)(0| + |0){l|) 



1 



E-/flead(l) + i(5 



(C17) 



Similar reasoning as the one described in Eqs. (|C5I 
(|C14p leads us to 



E-t^Gl^^{E'^)^t5 



(C18) 



with 



GLd(^;i) = (^ = i| 



^-i^lead(l)+^^ 



\x = l). (C19) 



Thanks to the translational invar iance, we should have 
GLd(^;0) = ^Ld(^; l)- Then, Eq. ((CHl) reduces to a 
quadratic equation 



which is followed by 



(C20) 



^Ld(^; 0) = ^ ^"If for -2t<E< 2t, 



2^2 



(C21) 



to^y E + iV4t2 - E 



(C23) 



Although we have derived the expression (|Cip partic- 
ularly for the present system (p!2|) with all the leads at- 
tached to a single site, the expression (|Cip itself holds for 
more general systems with appropriate changes of defi- 
nition of the effective Hamiltonian (|C2P : see Refs. [2I. IsHI. 
[m, [9^. We can reduce the calculation of the Green's 
function further for the present system (p!2|) , using the re- 
solvent expansion (|C13p again. For i = j = 0, Eq. (|C13p 
now gives 

+ • • • , (C24) 



where 



gUe) ^ {do\ 



'E-Hd 
Summing the series we obtain 



1^0). 



1 



Km 



(C25) 



(C26) 



This reduces the calculation of Gqq from inversion of an 
non-Hermitian matrix {E — H^) to a Hermitian matrix 
{E — Hd). For j = or i = 0, we have 



1 - ^ooKff ' 



/^d 



^d T/R' 



(C27) 
(C28) 
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respectively, and for general i and j we have 



' 1 r<d T/R ^Oj • 



(C29) 



Now we show how we can calculate all resonant states 
for the system ([12]). As is evident in the Fisher-Lee re- 
lation dH]), the conductance of the present system has 
poles in the complex energy plane wherever the Green's 
function G^q{E) has poles. The expression ()C26p imme- 
diately gives the equation for the resonant states in the 
form 



Ggo(^)Ki(^) = l. 



(C30) 



The resonant states given for the examples in Sec. IV 
were thus calculated. Equation (|C3Qp holds particularly 
for the present system (p!2]) with all the leads attached to 
a single site. For more general case, the Green's function 
Gqq is given by inversion of the matrix {E — H^{E)). 
Therefore, all resonant states can be calculated by solving 
the equation 



det{E - HfsiE)) = 0. 



(C31) 



The above discussion leads us to a much simpler way of 
deriving Eq. (jCip [99]; we formulate the Green's function 
so that it may have poles for the resonant states. Since 
a resonant state satisfies the boundary condition ([6]), we 
have 



(C32) 



with Re fc''®^ > 0. This terminates the Schrodinger equa- 
tion for the semi-infinite leads with the effective potential 



(C33) 



Solving the dispersion relation Ek = —2t cos k = — t(e* 



-ik 



), we have 



Ak 



2t 



(C34) 



which again gives Eq. (|C4p . See Ref. [99] for details. 



where G^q = (^01^^(^)14), Aqo = {do\A{E)\do), Too = 
{do\T{E)\do) and we made use of the fact that the ma- 
trix r has only the (0,0) element, Foo, for the present 
system (p!2]) : see Eq. (j32j). The solution is given by 



00 



Aqo ^ 1 
2 iFoo 




■ (D2) 



We here use Eqs. (|29|) and (|3T|) for z = j = 0, which gives 

(D3) 
(D4) 



Aoo — Gqo + G^Q, 



•-P _ -^00 ~ '-^OO 

00 - r^r^ ' 
'-^oo'-^oo 



The solution (|D2j) then may be written as 



^R 

'-^oo 



^00 y 



2 (^00 ~ G^o) 



(D5) 



This gives a consistent solution to Eq. ()D2p if we choose 
the ± sign as follows: 



Choose if (G^o)' + (G^o)' < 0; 
Choose if (G^o)^ + (Goo)^ > 0. 



Since Gqq = (G^q)* ^^r real E^ we have 
and hence we can write Eq. (jD6p as follows: 



Choose if |ReG^o| < llm^ail ; 
Choose "-" if |ReG^o| > |lmG^o| • 



(D6) 
(D7) 



(D8) 



(D9) 
(DIO) 



By using the expression (|C26p of the Green's function, 
we can further reduce the above criterion as follows: 



APPENDIX D: CHOOSING THE SIGN OF THE 
SOLUTION TO THE RICCATI EQUATION 

In this appendix we will derive a criterion to choose 
either the plus or the minus sign in the solution to the 
matrix Riccati equation, Eqs. (|38|) and (|39|) . We will 
consider the matrix element Gqq{E) = {do\G^ {E)\do) , 
which appears in the conductance (pij) . For i = j = 0, 
the equation ([36|) reduces to 



iG^orooG^o + G^o(2 + ^rooAoo) + Aqo = 0, (Dl) 



Choose if 



Choose "-" if 



(Ggo)-^-fE(T 

a 

(GS.)--fE(T 



J- 00 

2 ' 
(Dll) 

Too 
2 

(D12) 



This is the criterion to choose the sign in Eq. (jD2j). Using 
Eqs. (|C27p - (|C29p . we can show that this same criterion 
applies to Eqs. ([36|) and ([37|) . 
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